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(N ! Abstract 

Q_i' We investigate the occurrence of vibrational resonance in both classical and quantum mechanical 

<£ ' - . 

Morse oscillators driven by a biharmonic force. The biharmonic force consists of two forces of 

in ■ 

widely different frequencies lo and f2 with f2 3> uj. In the damped and biharmonically driven 
Q \ classical Morse oscillator applying a theoretical approach we obtain an analytical expression for 

the response amplitude at the low- frequency u. We identify the conditions on the parameters for 
the occurrence of the resonance. The system shows only one resonance and moreover at resonance 
the response amplitude is \/{duS) where d is the coefficient of linear damping. When the amplitude 

of the high-frequency force is varied after resonance the response amplitude does not decay to 

00 

zero but approaches a nonzero limiting value. We have observed that vibrational resonance occurs 

OV 

when the sinusoidal force is replaced by a square-wave force. We also report the occurrence of 



resonance and anti-resonance of transition probability of quantum mechanical Morse oscillator in 
the presence of the biharmonic external field. 
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I. INTRODUCTION 

Amplification and detection of weak signals are important in many branches of science. 
During the past two decades, a great deal of interest has been focused on signal processing in 
nonlinear systems. There are few interesting ways of enhancing the response of a nonlinear 
system to a weak signal. In stochastic resonance an optimum weak noise amplifies the 
response of a nonlinear system making use of the bistability of the system [T|2] , In another 
approach the noise is replaced by a high-frequency periodic signal. That is, the system is 
essentially driven by two periodic signals of widely different frequencies, say, u and Q with 
Q 3> u>. When the amplitude of the high-frequency external force is varied the oscillation 
amplitude of the system at the low-frequency u exhibits resonance. This high-frequency 
force induced resonance is termed as vibrational resonance [3}(5] . One-way or unidirectional 
coupling can also be able to improve the performance of coupled systems [6HH]. Similarly, 
the same phenomenon can be also observed when besides a noise or a high-frequency signal, 
a chaotic signal is used to perturb the system [9lfT0]. 

The study of vibrational resonance in the presence of a biharmonic force has received 
a considerable interest in recent years after the seminal paper by Landa and McClintock 
[3]. The occurrence of vibrational resonance has been analyzed in systems with monostable 
[TT] . bistable [2H5], spatially periodic states [12], excitable systems fT3][T4] . coupled oscillators 
[T5] and small world networks [T4|16|ll7j . Experimental evidence of vibrational resonance in 
vertical cavity laser system [T81I20] and in an electronic circuit [2T] has also been reported. 
Frequency-resonance-enhanced vibrational [22] and undamped signal propagation in one- 
way coupled oscillators [23] and maps [23] assisted by biharmonic force are found to occur 
as well. 

It is important to investigate the vibrational resonance in different kinds of systems and 
bring out the various features of it and the influence of characteristics of the potential of the 
system on it. In this connection we point out that (i) one or more resonances in monostable 
systems [TT] . (ii) a sequence of resonance peaks in a spatially periodic potential system [12] . 
(iii) additional resonances due to asymmetry in the potential [25], periodic and quasiperiodic 
occurrence of resonance peaks in time-delayed feedback systems [26H29] and large amplitude 
of vibrational resonance at the dynamic bifurcation point [30] have been reported. 

Motivated from some of the previous results, in the present paper we report our inves- 
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FIG. 1: The shape of the Morse potential for three values of f3. 



tigation on the vibrational resonance in the Morse oscillator. The potential of the Morse 
oscillator is 

V(x) = ~/?e-* (e- - 2) , (1) 

where /3 is a constant parameter representing the dissociation energy. Figure [I] depicts the 
form of the potential for a few values of j3. The potential is nonpolynomial and V(x) — y oo 
as x — y — oo while it becomes in the limit of x — > oo. It has one local minimum at x = 
and the depth of the potential is 0/2. The Morse oscillator was introduced as a useful model 
for the interatomic potential and fitting the vibrational spectra of diatomic molecules. It 
is also used to describe the photo-dissociation of molecules, multi-photon excitation of the 
diatomic molecules in a dense medium or in a gaseous cell under high-pressure and pumping 
of a local model of a polyatomic molecule by an infrared laser [3TH34] . The existence and 
bifurcations of periodic orbits have been studied in detail [5o"15o] . 

The equation of motion of the damped and biharmonically driven classical Morse oscil- 
lator is given by 



x 



+ dx + (3e x (l - e x ) — f cos ujt + g cos Qt, Q > u>. 



(2) 



Because of the difference in the time scales of the two periodic forces the motion of the 
system ()2]) contains both a slow component X(t) with period T = 2tt/u and a fast varying 



component i/j(t, r = Qt) with period 2n/Q. That is, we assume that x(t) = X(t) + ip(t, r). 
When the amplitude g or frequency Q is varied the response of the system fl2]) at the frequency 
u can exhibit resonance at a particular value of the control parameter g or Q. To analyse this 
resonance phenomenon and the influence of the shape of the potential on the resonance we 
obtain an equation of motion for the slow variable applying a perturbation theory. From the 
solution of the linearized version of the equation of motion about its equilibrium point, we 
find the analytical expression for the response amplitude Q which is the ratio of the amplitude 
of oscillation of the output of the system at the frequency u and the amplitude / of the input 
signal. From the expression of Q, we extract various features of vibrational resonance in the 
Morse oscillator. Particularly, we determine the value of g at which resonance occurs, the 
maximum value of Q at resonance and the limiting value of Q. We confirm all the theoretical 
predictions through numerical simulation. The theoretical treatment used for the analysis 
of the vibrational resonance with the periodic force / cos ut + g cos Qt can be applied for 
other types of biharmonic forces. In particular, we illustrate this for a square-wave form of 
low-frequency and high-frequency forces. Next, we consider the quantum mechanical Morse 
oscillator subjected to the biharmonic external field. Applying a perturbation theory we 
obtain an analytical expression for the first-order transition probability Pfi for a transition 
from an zth quantum state to an /th quantum state in time T caused by the applied external 
field. We analyse the influence of the high-frequency field on Pfi and show the occurrence 
of resonance and anti-resonance. 

II. CLASSICAL MORSE OSCILLATOR 

In this section we consider the classical Morse oscillator described by the equation of 
motion (J2J) and analyse the occurrence of vibrational resonance. 

A. Theoretical approach 

To find an approximate solution of Eq. (J2J, we make use of the method of separation of 
variables by assuming x = X(t) + ip(t, r — Qt), where X and ip are slow and fast variables 
respectively. Substituting x = X + ip in Eq. (J2J) and adding and subtracting the terms 



/3{e ^}e x and /3(e 2 ^)e 2x where (it) = (l/2ir) f Q ™ u dr, r = Vlt we obtain 

X + dX + /3e" x ((e"^) - (e^e - *) = / cos wt, (3) 

$ + # + f3e~ x (e~^ - (e^)) 

-{3e~ 2X (e~ 2 ^ -{e^}) = gcosM. (4) 

Because, ?/> is a rapidly changing variable with fast time r we approximate Eq. 03]) as ^ = 
gcosQt. The solution of this equation is ^ = ficosQt where /i = g/Q 2 . This solution gives 

<e-*> = ±J\» COST dT = I (f,), (5a) 

(e"*) = 7o(2/*), (5b) 

where Iq{z) is the zeroth-order modified Bessel function [57] . Now, Eq. ([3]) becomes 

X + dX + /3e~ x (J (/x) - J (2/i)e~ x ) = /coswt (6) 

Equation ([6]) can be treated as the equation of motion of a particle experiencing the external 
periodic force / cos ut and linear friction force in the effective potential 

V cS = ^f3e~ x [Wfi)e- X - 2/ (//)] . (7) 

In addition to the parameter fl, the effective potential depends on the parameters g and Q. 
The dependence of V e s on g and Q is in the form of a zeroth-order modified Bessel function. 
Figure [2(a) depicts the variation of io(/-0 an d Iq(2[x) with g for fi = 10. As g increases 
from zero, the values of io(/-0 an d Io{2fj) increase from the value 1. Io(2fi) increases more 
rapidly than 7 (yu) with g. Figure[2](b) shows the change in the shape of the effective potential 
with g. The location of the minimum of the potential V e g(X) or the X-component of the 
equilibrium point of the system ([6]) in absence of the periodic driving force is given by 

X* = -ln(44V (8) 

Because Io(2fi) > Io(fj), X* is always > and it moves away from the origin as g increases. 
Consequently, the depth AV c s =\ V e s(X*) | = /3/Q(/i)/2J (2/i) decreases with increase in 
g and the potential V e s becomes more and more flat. Later, we point out an important 
consequence of this on the response amplitude Q. 
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FIG. 2: (a) Dependence of a zeroth-order modified Bessel function with the control parameter g 
with SI = 10 and fi = g/O, 2 . (b) Plot of the effective potential V e s given by Eq. (J7J) for three values 
of g with 13 = 2 and U = 10. 

A slow oscillation takes place about X*. Therefore, for convenience, we introduce the 
change of variable Y = X — X* so that the slow oscillation occurs around Y* = 0. In terms 
of Y, Eq. 06]) becomes 



Y + dY + w r V y (1 - e~ y ) = / cos ut, 



where 



uj 2 = (3 



i 2 M 



(9a) 



(9b) 



For |/| ^C 1 it is reasonable to assume that the amplitude of Y is small so that we write 
series expansions for e~ Y and e~ 2Y and neglect the nonlinear terms in Y. This results in the 
linear equation 

Y + dY + u 2 T Y = f cos ut. (10) 

u r is the resonant frequency of oscillation of the slow motion. In the long time limit, the 
solution of Eq. ffTUj) is Y = Qf cos(ut + (p) where 
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s={u 2 -u 2 y + d 
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FIG. 3: (a) Response amplitude Q versus the control parameter g for the Morse oscillator for three 
values of /3. The values of the parameters are d = 0.5, / = 0.1, u = 1 and SI = 10. The values 
of j3 for the curves 1, 2 and 3 are 2, 1.5 and 0.5 respectively. The continuous curve and solid 
circles represent theoretically and numerically computed values of Q respectively. The horizontal 
dashed line denotes the limiting value of Q, Ql(<? — > oo). (b) Dependence of cj;? with g for the 
three values of /? used in the subplot (a). The solid circles on the g axis denote the values of g at 
which resonance occurs. 



where Q is the response amplitude of the system 
signal. 



at the low-frequency u of the input 



B. Analysis of the vibrational resonance 



In order to verify the theoretical predictions to be obtained from the analysis of Q given 
by Eq. (TTTK). we calculate Q from the numerical solution of Eq. ([2]). We use the formula 
Q = VQs + Qc/f where [31I5TT3] 

(12a) 

2 r T 
Q c = — — / x(t) cos u>t dt, (12b) 



Q s = — / x(t) smut dt, 

nT J Q 



with T = 2tt/u and n is big enough, say 500. 



The value of a control parameter at which resonance occurs (Q becomes maximum) 
corresponds to the minimum of the function S given by Eq. (ITTb ). The following are the 
key results of analysis of the theoretical expression of the response amplitude Q. 

• For a fixed value of g when u is varied resonance occurs when u = u VR given by 
(obtained from dS/duj = 0) 



d 2 d 2 

U VR = \JUt ~ J' U 2 > —. (13) 

When g is varied the condition for resonance is dS/dg = 4u r u rg (u 2 — u 2 ) = 0. Because 
Jo is always > and increases monotonically with g we have u r ^ and u; rg = duj T /dg ^ 
and hence the resonance condition is uj 2 = u 2 , that is, (3lQ(g/Q 2 )/I (2g/Q) = u 2 . 
Resonance occurs whenever the resonant frequency u r matches with the low-frequency 
u> of the periodic force. 

Figures [3]^a) and (b) show the variation of the response amplitude Q and u 2 respec- 
tively with g for three fixed values of (3 and d = 0.5, / = 0.1, u> = 1 and Q = 10. 
The theoretical value of Q is in very good agreement with Q obtained from the nu- 
merical solution of Eq. ([2]). In Fig. [3](a) for both /3 = 2 and 1.5 as g increases from 
zero the value of Q increases monotonically, it reaches a maximum at g = g VR and 
then decreases. For /3 = 2 theoretical and numerical values of g VR are 176 and 173 
respectively. In Fig.[3]^b) for both (3 = 2 and 1.5 at g = g VR (indicated by solid circles), 
u 2 = u 2 . 

At g = 0, Io(fi) = Iq (2jm) = 1 and u 2 (g = 0) = (3. As g increases io(/-0 an d Io(2fi) 
increases rapidly with Iq(2/j,) growing faster than Iq{^) as shown in Fig. [2](a). Conse- 
quently, u 2 decreases rapidly from the value of (3 for a while and then decays to zero 
slowly. The maximum value of u 2 is (3 and this happens at g = For (3 < u 2 , u 2 is 
always less than u 2 , that is, u 2 ^ u 2 for g > implying no resonance. This is shown 
in Fig. |3]^a) for /3 = 0.5 and u = 1 for which Q decreases continuously with g. 



Resonance is possible only for (3 > to . In this case, as g increases from zero, to 2 
decreases and becomes u 2 at only one value. Hence, there is only one resonance. This 
is shown in Fig. |3]^a) for (3 = 1.5 and 2 with u = 1. Figure H] presents the dependence 
of Q on (3 and g for to = 1. Resonance is seen only for (3 > u 2 {= 1). 
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FIG. 4: Three-dimensional plot of the theoretically computed response amplitude Q as a function 
of (3 and g. 
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FIG. 5: Variation of theoretically predicted (continuous curve) and numerically computed (solid 
circles) g VR with the parameter f3 for a few fixed values of the parameter u>. The value of 0, is 10w. 

• At resonance u 2 = u 2 and hence (5 max = l/(dui) and it depends only on d and uj. The 
response amplitude at resonance is independent of the parameters /3, g and Q. 

• An analytical expression for g VH (at which Q becomes maximum) is difficult to obtain 
because u 2 is a complicated function of g. However, g VR can be calculated from the 



resonance curve. It depends on fl, uj and (3 and independent of d and /. Q decreases 
with increase in the value of d. In Fig. [5] we plot the theoretical g VR and numerically 
computed g VR versus the parameter j3 for a few fixed values of uj. g VR increases with 
increase in the value of /3. 

• For very large values of g, uj 2 — > and Q approaches the limiting value Ql given by 

Q,(g ^oo)= ^== . (14) 

ujyuj 2 + d z 
That is, Q does not decay to zero but approaches the above limiting value (see 

Fig. [3(a)). This limiting value depends only on the parameters uj and d (note that 

Qmax = ~L/(duj)). The point is that when uj 2 — > 0, Eq. (J5J) becomes the damped free 

particle driven by the periodic force whose solution is 

Y(t) = QJcos(ut + $), $ = tan- 1 (rf/w). (15) 

For the Duffing and quintic oscillators [3|5|ll|25j V(x) (as well as the effective 
potential)— > oo as x — > ±oo. In this case the resonant frequency diverges after a 
few oscillations and thus Q decays to zero for large values of g. 

III. RESONANCE WITH A SQUARE-WAVE SIGNAL 

In this section we show that vibrational resonance can be realized in the Morse oscillator 
when the external signal is a square-wave and the theoretical analysis employed in the 
previous section for sinusoidal force can be applied to this case also. 

In the system (J5J) in place of iq = / cos ut + g cos Qt we consider three other forms of 
external force given by 

^(t) = f cos ut + g sgn(cos fit), (16) 

Fs(t) — f sgn(cos out) + g cos fit, (17) 

F 4 (t) = fsgn(cosut) + g sgn(cos fit) , (18) 

where sgn|w| denotes sign of u. In the theoretical analysis we use a Fourier series expansion 
for the square- wave signal. For the Morse oscillator driven by the periodic force Fi(t) the 
solution of slow motion is Yi = Q2/ cos(o;t + 02 ) where 

Q 2 = l ^=, (19a) 

J(ul 2 - UJ 2 ) 2 + d 2 UJ 2 

10 



2 - U(jjl) = — / e"^ 2 dr, (19b) 



where 

^ (-l)^cos(2n + l)r (igc) 

vr^ (2n+l)» • li9CJ 

For the case of the system driven by the force F3 the solution Y$ is 

Y, = *L± <=^ cos[(2, i + 1 M + fe ] r , (20.) 

n=0 



(a; r 2 3 - (2n + l) 2 w 2 ) + (2n + l) 2 d 2 



where 

" r ' 3 " ToM ' (20b) 

The solution Y3 has frequencies u and odd integer multiples of it. The response amplitude 
Qs corresponding to the fundamental frequency uj is 

4 1 

Qs = - , 9 =■ (21) 

VK3-W 2 ) 2 + ^ 

When the biharmonic force is chosen as _F 4 given by Eq. (fl8|) then the expression for the 
solution Y4 is the same as Y 3 except that now u 2 3 is replaced by a; 2 4 where 



2 PU 2 (V) 



<4 = TTT^V- (22) 

Then 



4 1 

Qa = = ■ (23) 



7T 



(w? >4 -o; 2 ) +d 2 c; 2 



Figure [6] shows both theoretically determined Q and numerically computed Q versus g 
with different input signals where d = 0.5, (3 = 2, f = 0.1, u; = 1 and Q = 10. In 
this figure, apart from close agreement of theoretical Q with numerical Q, we notice that 
Qi.max = Q2,max, Qs.max = Q4,max, Qi,l = <?2,l (limiting values of Q in the limit of g ->• 00) 
and Q 3j l = (?4,l- These results can be accounted from the theoretical expressions of Qj's. 
The resonance condition u 2 { = u 2 , i = 1, 2, 3, 4 in the expression of Qj's gives 

*Vl,max *V2,max ~5 j I ^48.) 

do; 
4 4 

*t;3,max V4,max J ^;l,max; I Z4DI 

TTdUJ IT 
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FIG. 6: Variation of the response amplitude Q with g for four different types of biharmonic forces. 
The continuous curve and the solid circles represent theoretically calculated Q and numerically 
calculated Q respectively. The biharmonic forces corresponding to the curves and solid circles 
marked as 1, 2, 3 and 4 are / cos cot + g cos Qt, F2, F3 and F4 respectively. 



that is, Q3 )max and Q4 imax ar e 4/7r ps 1.27324 times of Qi <m£LX . For the parametric values 
used in our analysis, Qi imax = 2 and hence Q3 jmax and Q^ mabX are 2.54648 as is the case in 
Fig. [61 For sufficiently large values of g, we have w^ ~ and hence 



1 4 

<2l,L = Q2,L = , = = , <53,L = Q4,L = — Ql,L- 

WVU 2 + rf 2 7T 



(25) 



Furthermore, because of w 2 x = a; 2 3 and w 2 2 = w 2 4 we find that 



_ 4 _ 4 

Q3 = — Qi 5 Q4 = — Q2- 

TX TV 



(26) 



That is, the response amplitude at the frequency u when the input signal is a square-wave 
with fundamental frequency being u is 4/ V times that of the signal / cos cut. Numerical 
results in Fig. [6] confirms all the above theoretical predictions. 

Our analysis shows that the form of the low- and high-frequency forces need not be 
identical. For any arbitrary force containing a component with frequency u, enhancement 
of the amplitude of the output signal at the frequency to can be achieved by using another 

12 



arbitrary force containing a frequency Q ^> u>. When the force involved are simple periodic 
function of t then the theoretical analysis of vibrational resonance is very much feasible. 

IV. QUANTUM MECHANICAL MORSE OSCILLATOR 

In the previous two sections we focused our analysis on vibrational resonance in the clas- 
sical Morse oscillator. In the present section we are concerned with the quantum mechanical 
Morse oscillator in the presence of the biharmonic external field W(t) = F cos ut + gcosQt 
with Q 3> u>. When a quantum mechanical system is subjected to a time-dependent exter- 
nal field the system undergoes transition between the energy eigenstates. Therefore, we are 
interested in knowing the probability of finding the system in an fth state at time t. 

The unperturbed Hamiltonian of the system is H = p 2 /(2m) + V(x) where V(x) is given 
by Eq. (CQ). The unperturbed system H Q (j) n = E n cj) n is exactly solvable and the eigenf unctions 
and energy eigenvalues are given by 



where 



<P n = N n z x - n - 1 ' 2 e- z l 2 L k n {z), (27a) 



and 



* = «•-. * = % M(^WF' (27b) 

k = 2A - 2n - 1, Ihz) = ^-^-e~ z z n+k (27c) 

nV ' n\ dz n v ' 

h 2 ( i\ 2 l 

E n = -—\\-n--\ , n = 0,l,--- andn< A--. (28) 

In Eq. (127b) L^(z) are the generalized Laguerre polynomials. The Morse oscillator has a 
finite number of bound states and the number of bound states can be controlled by the 
parameter (3. Setting the values of ft and m as unity for convenience and /3 = 9 we obtain 

E = -3.125, E 1 = -1.125, E 2 = -0.125. (29) 

There are only three bound states. Figure [7] shows the energy eigenvalues and the eigen- 
functions for (3 = 9. 

In the presence of the external field the Hamiltonian of the system is H = Hq + eHi, 
where Hi = xW(t). The time-dependent Schrodinger equation for the system is 

ih^- = Hi/,, (30) 
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FIG. 7: Bound state energy eigenvalues and eigenfunctions of the Morse oscillator for f3 = 9. The 
dashed lines are the energy eigenvalues. The potential V(x) is also shown. 



where ip(x,t) is the wave function of the perturbed system. We write 



ip(x,t) = ^a n (t)0 n (x)e ' 



iE n t/h 



(31) 



The probability of finding the system in the state n is P n (t) = |a n (t)| 2 , ^ n |a n (t)| 2 = 1. 

To determine a n (t) we apply the standard time-dependent perturbation theory pET] . Sup- 
pose the external field is switched-on at t = and switched-off at t — T, that is, the external 
field is applied during a finite time interval T. Assume that the system is initially in a state 
% with the eigenfunction fa. Then at t = the probability of finding the system in the state 
i is 1 and the probability of finding the system in the other states is 0: o n (0) = 5 n i. Due to 
the applied field the system can make a transition from the state % to an another state after 
the time T. Once the perturbation is switched-off the system settles down to a stationary 
state and denote this final state as /. 

To determine a n (t) we substitute Eq. ( |3T1) . aj = a\ + ea\ + • • • and equate the terms 
containing various powers of e to 0. Up to first-order in e, after some algebra, we obtain 



a f = Sfi and 



.W 



C 



/'■ 



af(T) = ^s 



(32a) 



14 




FIG. 8: log |s| 2 versus w for the states / = 0, 1 and 2 for j3 = 9, F = 0.05, # = 0, T = tt and i = 0. 



where 



n± 



Wfi 



F(r 1+ + n_) + g(r 2+ + r 2 _), 

1 _ gi(w /i ±a;)T ]_ _ e i(w /8 ±C)T 

?"2± = 



W/j ± w 



Ufi ± n 



/oo 
0}x0j dx. 
-oo 



(32b) 
(32c) 

(32d) 



,(i) ( 



The transition probability for zth state to /th state is given by PfiiT) = 5fi + ea\ (T) 
In a, (T) the term s alone depends on the parameters F, u>, g and f2 of the external field 
and T. Therefore, we study the variation of the quantity |s| 2 with the parameters of the 
external field. 

We fix /3 = 9, F = 0.05, T = it and assume that the system is initially in the ground state 
(i = 0). Figure [8] shows the variation of log |s| 2 with u for g = 0. The first-order correction 
to transition probability displays a sequence of resonance peaks with decreasing amplitude. 
The following results are evident from Eqs. f)32b-c) and Fig. |HJ The quantity s consists of 
only r i+ and ri_. We have cuoo = 0, uio = 2 and W20 = 3. Consequently, for the states 
/ = 0, 1 and 2 the quantity r i+ can be neglected when uj ~ 0, 2 and 3 respectively because 
the denominator in ri_ is ~ 0. Thus, the first-order transition probability for the states 
/ = 0, 1 and 2 becomes maximum at u> = 0, 2 and 3 respectively. For u> rs the values of 
|s| 2 for / = 0, 1 and 2 are ~ 4F 2 ti 2 , and 16F 2 /9 respectively. The first-order transition 
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probability for the state / = 1 is ~ 0. For the / = state |s| 2 = when u) = 1, 2, • • • and it 
becomes maximum when u — n + |, n — 1, 2, • • • . In the case of / = 1 at u = 4, 6, ■ ■ • the 
numerator in both r± + and ri_ becomes zero and thus the quantity |s| 2 is minimum. For 
/ = 2, at odd integer values of a; except at u = 3, r±± = and hence \s\ 2 becomes minimum. 

Next, we include the high-frequency field and vary its amplitude g and the frequency fl. 
Figure M presents the results for a few fixed values of u> with Q = 5w. When oj = 1, in 
Eq. fl32b). r 1+ + ri_ = and r 2+ + r 2 _ = for / = and 2 and hence the increase of g 
has no effect on |s| 2 . For the / = 1 state we find s = (8/3) .F — (8/21)g. As g increases 
|s| 2 decreases from (8/3)F, becomes at g = 7F(= 0.35) and then increases with further 
increase in g as shown in Fig. E(a). The first-order transition probability of the state / = 1 
exhibits anti-resonance at g = 0.35. Anti-resonance can be realizable for other states also 
for appropriate choices of u. For example, in Figs. E^b) and (c) corresponding to u — 1.7 
and 2 respectively we can clearly notice anti-resonance for / = state and / = 2 state. In 
Fig. [9jd) where u = 3.5 the quantity \s\ 2 increases monotonically with the control parameter 
g for all the three states. 

Because r 2 + and r 2 _ (given by Eq. ()32b)) contain terms which are sinusoidal functions 
of Q the first-order transition probability can exhibit a sequence of resonance peaks when 
Q is varied for fixed values of the other parameters. This is shown in Fig. [10] for four sets 
of values of u and g. In all the cases Q is varied from 5a;. In Figs. [TOT a-c) |s| 2 = at the 
starting value of Q for the / = state. |s| 2 ^ for wide ranges of values of Q. Transition 
probability of all the states show a sequence of resonance peaks. In Fig. [TOT d) \s\ 2 of / = 
state is close to its values of the other two states. However, as Q increases |s| 2 of the / = 1 
and / = 2 states oscillate but ^ 0. In contrast to this, |s| 2 of the / = state becomes at 
certain values of Q. 

It is noteworthy to compare the effect of high-frequency external force in the classical and 
quantum mechanical Morse oscillators. A classical nonlinear system can exhibit a variety 
of dynamics when a control parameter is varied. However, in the Morse oscillator for the 
choice |/| <C 1 and Q ^> u> when the amplitude g of the high-frequency force is varied 
the system is found to show only periodic motion with period T = 2ix/oj. The response 
amplitude of the motion exhibits a single resonance when the control parameter g or u or 
Q is varied. Resonance occurs whenever the resonant frequency u r (Eq. ([9b)) matches with 
the frequency u. In the case of the quantum mechanical Morse oscillator we have considered 
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FIG. 9: Variation of \s\ 2 (in log scale) with the amplitude g of the high-frequency external field for 
four fixed values of the frequency 00 of the low- frequency external field. The value of £1 is fixed as 
huj. 

the simple case of switching-on the external field at t — and switching-off it at t — T. 
In the absence of the high-frequency field the first-order transition probability Pfi shows a 
sequence of resonances with decreasing amplitude when the parameter u is increased. The 
dominant resonance occurs at u — Ufi. Resonance is not observed when the amplitude g 
of the high-frequency field is varied. However, anti-resonance of Pfi takes place for certain 
values of u. Multiple resonance of Pfi occurs when the frequency Q of the high-frequency 
field is varied. 

V. CONCLUSION 

We have reported our investigation on high-frequency periodic force induced resonance 
at the low-frequency component of the output of a single oscillator. Using a perturbation 
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FIG. 10: log |s| 2 versus $7 (the frequency of the high-frequency external field) for four sets of values 
of u> and g. 



theory an analytical expression for the response amplitude is obtained. Interestingly, from 
the analytical expression of Q we are able to derive various features of the vibrational 
resonance and its mechanisms. The occurrence of the resonance depends on the parameter 
(3 and uj while the values of the response amplitude at resonance and for large values of 
g depend only on the damping coefficient d and the low-frequency u. In the system the 
limiting values of Q is nonzero because for sufficiently large values of g, u^ ~ 0. From the 
analytical expression of Q given by Eq. lTTTk ) we note that Q becomes a nonzero constant if 
bJl ~ a constant. The theory used in our present analysis can be applied to the other different 
forms of input signal than coswt and sin cut. We have demonstrated its applicability for the 
case of square-wave form. All the theoretical results are well supported by the numerical 
simulation. We have considered the quantum version of the Morse oscillator in the presence 
of the biharmonic external field. Interestingly, a high-frequency external field is found to 
induce resonance and anti-resonance on the transition probability for the transition from an 
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ith state to an /th state. If the transition probability P/j of a state is weak in the presence 
of a harmonic external field with a particular frequency u, then it can be enhanced by an 
another external field of relatively high-frequency. The dominance of a state can be changed 
by the high-frequency external field. That is, Pfi can be controlled by a second harmonic 
external field. 
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